Journal  of  Clinical  Microbiology,  Aug.  2011,  p.  2786-2797 
0095-1137/1 1/S  12.00  doi:  10. 1128/JCM.00033-1 1 

Copyright  ©  2011,  American  Society  for  Microbiology.  All  Rights  Reserved. 


Vol.  49,  No.  8 


Francisella  tularensis  Molecular  Typing  Using  Differential 
Insertion  Sequence  Amplification7 

Marilynn  A.  Larson,1*  Paul  D.  Fey,1  Amanda  M.  Bartling,1  Peter  C.  Iwen,1  Michael  P.  Dempsey,2! 

Stephen  C.  Francesconi,3  and  Steven  H.  Hinrichs1 

Department  of  Pathology  and  Microbiology,  University  of  Nebraska  Medical  Center,  Omaha,  Nebraska Division  of  Microbiology , 
Armed  Forces  Institute  of  Pathology,  Washington,  DC2;  and  Naval  Medical  Research  Center,  Silver  Spring,  Maryland 3 

Received  5  January  2011/Returned  for  modification  29  March  2011/Accepted  16  May  2011 

Tularemia  is  a  potentially  fatal  disease  that  is  caused  by  the  highly  infectious  and  zoonotic  pathogen 
Francisella  tularensis.  Despite  the  monomorphic  nature  of  sequenced  F.  tularensis  genomes,  there  is  a  significant 
degree  of  plasticity  in  the  organization  of  genetic  elements.  The  observed  variability  in  these  genomes  is  due 
primarily  to  the  transposition  of  direct  repeats  and  insertion  sequence  (IS)  elements.  Since  current  methods 
used  to  genotype  F.  tularensis  are  time-consuming  and  require  extensive  laboratory  resources,  IS  elements  were 
investigated  as  a  means  to  subtype  this  organism.  The  unique  spatial  location  of  specific  IS  elements  provided 
the  basis  for  the  development  of  a  differential  IS  amplification  (DISA)  assay  to  detect  and  distinguish  the  more 
virulent  F.  tularensis  subsp.  tularensis  (subtypes  A.I  and  A.II)  and  subsp.  holarctica  (type  B)  strains  from  F. 
tularensis  subsp.  novicida  and  other  near  neighbors,  including  Francisella  philomiragia  and  Franciselia-iike 
endosymbionts  found  in  ticks.  Amplicon  sizes  and  sequences  derived  from  DISA  showed  heterogeneity  within 
members  of  the  subtype  A.I  and  A.II  isolates  but  not  the  type  B  strains.  These  differences  were  due  to  a  312-bp 
fragment  derived  from  the  IS  element  ISFtul.  Analysis  of  wild-type  F.  tularensis  isolates  by  DISA  correlated 
with  pulsed-field  gel  electrophoresis  genotyping  utilizing  two  different  restriction  endonucleases  and  provided 
rapid  results  with  minimal  sample  processing.  The  applicability  of  this  molecular  typing  assay  for  environ¬ 
mental  studies  was  demonstrated  by  the  accurate  identification  and  differentiation  of  tick-borne  F.  tularensis. 

The  described  approach  to  IS  targeting  and  amplification  provides  new  capability  for  epidemiological  inves¬ 
tigations  and  characterizations  of  tularemia  source  outbreaks. 


Francisella  tularensis  is  a  facultative  intracellular  bacterium 
and  the  causative  agent  of  the  zoonotic  disease  tularemia  (10). 
This  Gram-negative  microbe  is  highly  infectious,  with  as  few  as 
10  organisms  being  capable  of  causing  disease  in  numerous 
mammalian  species,  including  humans.  Accordingly,  F.  tular¬ 
ensis  is  considered  a  potential  biological  weapon  and  is  classi¬ 
fied  by  the  Centers  for  Disease  Control  and  Prevention  (CDC) 
as  a  select-agent  pathogen  (29).  The  three  widely  accepted 
subspecies  of  F.  tularensis  include  F.  tularensis  subsp.  tularensis 
(also  known  as  type  A),  F.  tularensis  subsp.  holarctica  (also 
known  as  type  B),  and  F.  tularensis  subsp.  mediasiatica, 
whereas  classification  of  the  fourth  subspecies,  F.  tularensis 
subsp.  novicida,  remains  in  contention  (13,  14,  20).  Differen¬ 
tiation  among  these  subspecies  is  important  due  to  their  vari¬ 
ation  in  virulence  and  for  the  tracking  of  tularemia  outbreaks. 
Although  F.  tularensis  subsp.  mediasiatica  was  reported  to  ex¬ 
hibit  a  moderate  degree  of  virulence  in  mammals,  the  geo¬ 
graphic  distribution  in  Central  Asia  and  limited  accessibility  to 
isolates  have  restricted  the  knowledge  of  this  subspecies’ 
pathogenicity  (4).  F.  tularensis  subsp.  novicida  is  generally  con¬ 
sidered  an  opportunistic  microbe,  whereas  F.  tularensis  subsp. 
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tularensis  and  F.  tularensis  subsp.  holarctica  are  capable  of 
causing  a  severe  form  of  tularemia.  Moreover,  type  A  strains 
(F.  tularensis  subsp.  tularensis)  can  potentially  cause  a  life- 
threatening  illness.  The  virulent  type  A  and  type  B  clades  show 
remarkable  geographic  associations:  F.  tularensis  subsp.  holar¬ 
ctica  is  prevalent  throughout  the  Northern  Hemisphere,  while 
F.  tularensis  subsp.  tularensis  is  restricted  primarily  to  North 
America  (25).  Recently,  subtypes  A.I  and  A.II  of  F.  tularensis 
subsp.  tularensis  have  been  correlated  with  certain  regions  in 
North  America  and  differences  in  virulence,  with  subtype  A.I 
infection  being  reported  to  be  more  often  fatal  (18).  There¬ 
fore,  the  rapid  detection  and  typing  of  F.  tularensis  are  critical 
for  epidemiological  tracking,  clinical  disease  management,  and 
biodefense  monitoring. 

Numerous  molecular  studies  have  been  evaluated  for  the 
detection  and  identification  of  F.  tularensis  from  biological 
reservoirs  such  as  lagomorphs  and  ticks.  However,  amplifica¬ 
tion  of  the  16S  rRNA  gene  and  other  targets  has  produced 
misleading  results  due  to  their  cross-reactivity  with  near  neigh¬ 
bors  and  Francisella- like  endosymbionts  known  to  exist  in  ticks 
(3, 11, 17).  In  general,  examination  of  the  population  structure 
of  F.  tularensis  at  the  genetic  level  has  required  more  robust 
and  time-consuming  approaches,  such  as  pulsed-field  gel  elec¬ 
trophoresis  (PFGE),  multilocus  variable-number  tandem-re¬ 
peat  analysis,  and  whole-genome  detection  of  single-nucle¬ 
otide  polymorphisms  (SNP).  The  value  of  these  studies  has 
been  significant.  For  example,  phylogenetic  analyses  of  repeats 
and  SNP  content  suggested  that  F.  tularensis  subsp.  holarctica 
originated  from  North  America  and  was  introduced  multiple 
times  into  Eurasia  (15, 33).  Furthermore,  comparative  genome 
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hybridization  and  sequencing  were  used  to  demonstrate  a 
unique  deletion  present  within  an  F.  tularensis  isolate  found  in 
an  Iberian  Peninsula  tularemia  outbreak  (2,  7,  8). 

Other  genetic  markers  with  the  potential  for  use  in  molec¬ 
ular  typing  have  been  identified.  Seventeen  contiguous  regions 
were  previously  identified  in  F.  tularensis  subsp.  holarctica  that 
were  dispersed  in  F.  tularensis  subsp.  tularensis  using  paired- 
end  sequencing  (8).  A  recent  genomic  comparison  between 
four  F.  tularensis  subsp.  tularensis  isolates  (SCHU  S4,  FSC  198, 
NE-061598,  and  FSC  033)  confirmed  that  these  subtype  A.I 
strains  are  clonal,  with  a  high  degree  of  chromosomal  synteny 
(23).  Nevertheless,  distinctions  among  these  strains  could  be 
made  based  on  inversions  and  rearrangements  that  were  pre¬ 
dictably  mediated  by  direct  repeats  and  insertion  sequence  (IS) 
elements  (23).  Additional  computational  analysis  suggested 
that  several  of  the  contiguous  regions  could  be  used  to  differ¬ 
entiate  between  type  A  and  type  B  strains.  However,  since  few 
wild-type  isolates  have  been  sequenced,  the  application  of  this 
approach  to  characterize  new  foci  of  tularemia  outbreaks  was 
unknown.  To  address  this  challenge,  a  molecular  assay  was 
developed  using  differential  insertion  sequence  amplification 
(DISA)  to  characterize  wild-type  F.  tularensis  isolates,  and 
results  were  compared  to  the  typing  patterns  obtained  with 
PFGE  using  two  different  restriction  endonucleases. 

MATERIALS  AND  METHODS 

Bacterial  strains.  Eighty-eight  wild-type  and  four  reference  strains  of  F.  tula¬ 
rensis  were  evaluated,  as  shown  in  Fig.  2  and  3.  The  wild-type  isolates  were 
obtained  from  multiple  sources,  including  public  health  and  diagnostic  labora¬ 
tories  in  the  United  States.  The  identities  of  all  F.  tularensis  isolates  were 
confirmed  according  to  the  testing  protocols  utilized  by  Laboratory  Response 
Network  reference  laboratories.  The  reference  /•'.  tularensis  strains  included 
SCHU  S4  (subtype  A.I),  WY96-3418  (subtype  All),  the  live  vaccine  strain 
(LVS)  (type  B),  and  F.  tularensis  subsp.  nnvicida  (Table  1).  Twenty-three  non- 
Fmncisella  species  representing  21  different  species  as  well  as  one  Franciselta 
pluloinimgia  strain  were  included  in  an  exclusivity  panel  to  assess  specificity 
(Table  1).  Select-agent  strains  were  transferred  to  the  University  of  Nebraska 
Medical  Center,  Omaha,  NE.  according  to  the  requirements  of  the  Select  Agent 
Program  as  outlined  in  the  Animal  and  Plant  Health  Inspection  Service/CDC 
form  2.  guidance  document  for  request  to  transfer  select  agents  and  toxins  (5). 
The  manipulation  of  viable  culture  material  was  performed  by  authorized  indi¬ 
viduals  within  a  biosafety  level  3  laboratory  certified  for  select-agent  work  by  the 
U.S.  Department  of  Health  and  Human  Services,  using  laboratory  biosafety 
criteria  as  described  by  the  National  Institutes  of  Health  and  the  CDC  (9).  All 
Francisella  isolates  were  grown  on  chocolate  agar  plates  (Kernel,  Lenexa,  KS) 
and  incubated  at  37°C  with  5%  C02  for  3  days  before  processing. 

DNA  isolation.  The  extraction  of  chromosomal  DNA  from  all  bacterial  species 
was  performed  by  using  CI'AB  (cetyltrimcthyl  ammonium  bromide)  according  to 
standard  procedures  (34).  For  the  environmental  samples,  ticks  and  insects  were 
homogenized  in  cold  phosphate-buffered  saline  (PBS)  within  containment  by 
using  a  Mini  Beadbeater-8  instrument  (Bio  Spec  Products,  Inc..  Bartlesville,  OK) 
according  to  the  manufacturer’s  recommendations.  DNA  was  then  isolated  from 
the  arthropod  and  insect  extracts  by  using  the  QLAamp  DNA  blood  minikit 
(Qiagen,  Valencia,  CA).  Pond  water  and  soil  samples  were  subjected  to  DNA 
extraction,  as  was  performed  for  the  tick  extracts  but  without  homogenization. 
The  DNA  from  all  environmental  samples  was  concentrated  following  extrac¬ 
tion,  and  approximately  1  ng  of  DNA  was  utilized  for  the  PCR  assays.  If  an 
inhibitor  was  detected,  the  DNA  was  diluted,  and  0.1  ng  of  DNA  was  used 
instead.  Demtacentor  and  Ambtyonuna  ticks  were  collected  from  canine  or  bovine 
sources  as  well  as  by  flagging  or  dragging  and  were  identified  to  species  level 
according  to  previously  reported  keys  (16,  35). 

PFGE.  Agarose-embedded  chromosomal  DNA  for  PFGE  was  prepared  and 
digested  with  the  restriction  endonucleases  Pmel  and  BamHI  (Feimentas,  Inc., 
Glen  Bumie,  MD)  as  previously  described  (12).  Agarose  gels  (1%,  wt/vol)  con¬ 
tained  Francisella  DNA  digested  with  Pmel  or  BamHI  and  were  run  for  17.5  h 
or  163  h,  respectively,  with  an  initial  switch  time  of  1.8  s  and  a  final  switch  time 
of  10.7  s.  Migration  profiles  of  the  restriction  fragments  were  normalized  to 
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TABLE  1.  Reference  strains  of  Francisella  and  non-Francisella 
strains  used  in  this  study 


Species  Strain  designation(s)r 


Francisella  tularensis  subsp.  tularensis “ . SCHU  S4,  NR-643* 

Francisella  tularensis  subsp.  tularensis " . WY96-3418,  NR -644* 

Franciselta  tularensis  subsp.  holarctica . ATCC  29684,  LVS, 

NR-14* 

Francisella  tularensis  subsp.  novicida “ . ATCC  15482,  Utah 

112,  NR-13* 

Francisella  pltiloiniragia . ATCC  25015 

Acinetobacter  baumannii . Wild  type” 

Bacillus  airopliaeuS1 . ATCC  9372 

Bntcella  canis . Wild  typec 

Burkholderia  cepacia . . . ATCC  25608 

Corynebaclcrittm  pseudodiplillieriticum . ATCC  10701 

Enterobacter  cloacae . ATCC  13047 

Enterococcus  faecalis . ATCC  29212 

Escherichia  coli . ATCC  35218 

Escherichia  coli . Wild  type' 

Haemophilus  injluenzae . ATCC  10211 

Klebsiella  pneumoniae . ATCC  BAA-1705 

Micrococcus  htteus . ATCC  49732 

Neisseria  lactamica . ATCC  23970 

Proteus  vulgaris . ATCC  49132 

Pseudomonas  aeruginosa . ATCC  10145 

Serratia  marcescens . Wild  type' 

Staphylococcus  aureus . ATCC  25923 

Staphylococcus  aureus . Wild  typec 

Staphylococcus  epidermirlis . ATCC  14990 

Streptococcus  pneumoniae . ATCC  6305 

Streptococcus  pyogenes . ATCC  12384 

Yersinia  pestiS1 . ATCC  19428,  NR-18* 

Wolbachia  persica . ATCC  VR-331 


"  Received  through  the  submission  of  APH1S/CDC  form  2,  request  to  transfer 
select  agents  and  toxins. 

*  Strain  designation  from  the  NIH  Biodcfcnsc  and  Emerging  Infections  Re¬ 
search  Resources  Repository  (BEI  Resources). 

'  Wild-type  strain  from  the  University  of  Nebraska  Medical  Center  stock 
collection. 

''  Previously  known  at  Bacillus  subtilis  var.  niger. 

r  Abbreviations:  LVS.  live  vaccine  strain;  ATCC,  American  Type  Culture 
Collection. 


Smal-digested  Staphylococcus  aureus  NCTC  8325  by  using  Bionumerics  software 
(Applied  Maths,  Inc.,  Austin,  TX).  Cluster  analysis  was  performed  by  using  the 
Dice  correlation  coefficient  and  the  unweighted-pair  group  mathematical  aver¬ 
age  (UPGMA)  clustering  algorithm. 

PCR  amplifications.  The  continuous  region  10  (CR10)  primer  set  consisted  of 
the  oligonucleotides  5'-GCTGATGATTCTAGCCTTAAAGAAG-3’  (CR10  C), 
5'-GCCAATCCTACA’nYATAGAACCTG-3'  (CR10  S),  and  5'-GCCTGGCA 
TAATTACTG 1 1  f  l  AGC-3’  (CR10  L).  The  CR16  primer  set  was  comprised  of 
the  primers  5’-CTTGTCAGAGTTGGAGTGAAGC-3'  (CR16  C),  5’-CCACT 
ACCTCGAATCTTACACAAAG-3'  (CR16  S),  and  5'-GTCAAATCACTTTC 
CTCTCGTTC-3'  (CR16  L).  PCR  amplifications  were  performed  by  using  the 
indicated  CR  primers  and  Platinum  DNA  polymerase  (Invitrogen,  Carlsbad, 
CA)  as  recommended  by  the  manufacturer.  The  cycling  parameters  were  1  cycle 
for  2.5  min  at  95°C  and  40  cycles  for  30  s  at  95°C,  1  min  at  60“G  and  3  min  at 
72°C.  followed  by  1  cycle  for  5  min  at  72°C.  PCR  products  were  resolved  in  a 
0.7%  (wt/vol)  agarose  gel,  stained  with  cthidium  bromide,  and  visualized  with 
LTV  light.  Representative  amplicons  were  cloned  into  pCR4-TOPO  (Invitrogen, 
Carlsbad,  CA)  and  sequenced  by  using  M13  forward  and  reverse  primers  by  the 
University  of  Nebraska  Medical  Center  DNA  Sequencing  Facility. 

RESULTS 

PCR  and  chromosomal  mapping  of  amplicons  derived  from 
control  strains.  Computational  analysis  of  the  different  con¬ 
tinuous  regions  (CRs)  previously  described  (8)  suggested  that 
two  loci,  namely,  CR10  and  CR16,  could  be  used  to  differen- 
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tiate  between  F.  tularensis  type  A  and  type  B  strains.  The  CR10 
and  CR16  primers  were  based  on  chromosomal  content  within 
these  two  loci  that  was  common  (C)  to  both  SCHU  S4  and 
LVS  or  specific  to  either  SCHU  S4  (S)  or  LVS  (L).  The 
location  of  the  CR  primers  and  the  predicted  amplicon  size 
were  initially  determined  by  in  silico  analysis  for  all  Francisella 
species  and  subspecies  in  which  the  genome  sequence  was 
available,  including  F.  tularensis  subsp.  mediasiatica.  PCR  eval¬ 
uations  for  each  oligonucleotide  pair  in  the  CR10  and  CR16 
primer  sets  (C,  S,  and  L)  were  then  performed  by  using  three 
well-characterized  and  sequenced  strains  of  F.  tularensis ,  spe¬ 
cifically  SCHU  S4  (GenBank  accession  no.  AJ749949;  subtype 
A.I),  WY96-3418  (here  referred  to  as  WY96)  (GenBank  ac¬ 
cession  no.  CP000608;  subtype  A.II),  and  LVS  (GenBank  ac¬ 
cession  no.  AM233362;  type  B). 

Preliminary  PCR  assessments  revealed  that  primer  pair 
CR10  C  and  CR10  S  and  primer  pair  CR16  C  and  CR16  S 
were  specific  for  SCHU  S4,  generating  1.2-kb  and  2.3-kb  PCR 
products,  respectively.  Conversely,  primer  pair  CR10  C  and 
CR10  L  and  primer  pair  CR16  C  and  CR16  L  amplified  a 
single  PCR  product  that  differed  in  size  for  WY96  (1.3  kb  and 
1.7  kb,  respectively)  and  LVS  (2.1  kb  and  0.8  kb,  respectively). 
The  amplicon  sizes  obtained  with  the  CR10  and  CR16  primers 
for  these  three  reference  strains  correlated  with  the  expected 
lengths  based  on  the  nucleotide  sequences  deposited  in  the 
GenBank  database  (Fig.  1).  These  results  confirmed  the  in 
silico  predictions  for  all  F.  tularensis  subspecies,  except  F.  tu¬ 
larensis  subsp.  mediasiatica ,  for  which  no  DNA  was  available. 

The  PCR  products  obtained  for  SCHU  S4  with  primer  pair 
CR10  C  and  CR10  S  and  for  WY96  and  LVS  with  primer  pair 
CR10  C  and  CR10  L  contained  the  IS  element  ISFtul  (Fig. 

IA) .  However,  an  approximately  300-bp  region  of  the  ISFtul 
gene  in  SCHU  S4  was  apparently  deleted,  unlike  the  presum¬ 
ably  intact  ISFtul  in  WY96  and  LVS.  The  absence  of  deoD  in 
the  SCHU  S4  CRIO-derived  amplicon,  the  truncated  ISFtul 
gene,  and  the  nucleotide  differences  in  primers  CR10  S  and 
CR10  L  indicated  that  a  unique  SCHU  S4  chromosomal  region 
was  amplified  compared  to  the  PCR  products  obtained  with 
WY96  and  LVS.  Furthermore,  an  additional  IS  element,  spe¬ 
cifically,  IS Ftu2,  was  present  in  the  amplicon  derived  from  LVS 
with  CR10  C  and  CR10  L  but  absent  in  the  related  WY96  PCR 
product.  These  findings  verify  that  this  chromosomal  region 
differs  between  subtype  A.II  and  type  B  strains  even  when  the 
same  CR10  primers  are  utilized  (Fig.  1A). 

For  the  amplicons  derived  from  primer  pairs  CR16  C  and 
CR16  S  as  well  as  CR16  C  and  CR16  L,  ISFtul  was  present  in 
SCHU  S4  and  WY96,  respectively,  but  absent  in  LVS  (Fig. 

IB) .  The  smaller  PCR  product  obtained  for  LVS  with  primers 
CR16  C  and  CR16  L  reflects  the  absence  of  this  IS  element. 
When  the  CR16-derived  chromosomal  regions  in  SCHU  S4 
and  WY96  were  compared,  ISFtul  and  the  adjacent  fnuDl 
gene  were  inverted,  with  an  apparent  duplication  of  fnuDl  in 
the  latter  strain.  Since  LVS  is  an  attenuated  strain,  the  chro¬ 
mosomal  contents  of  the  amplified  CR10  and  CR16  regions 
were  examined  in  all  available  F.  tularensis  subsp.  holarctica 
genomes.  These  assessments  indicated  that  the  targeted  re¬ 
gions  shared  considerable  nucleotide  sequence  identity  in  type 
B  isolates.  Collectively,  the  chromosomal  mapping  of  the 
CR10  and  CR16  primers  in  the  three  control  strains  concurs 
with  the  amplicon  lengths  obtained  and  confirms  that  the  am¬ 


plified  regions  differ  in  contents  among  the  F.  tularensis  sub- 
type  A.I,  subtype  A.II,  and  type  B  strains. 

PFGE  genotyping.  PFGE  analysis  was  shown  previously  to 
type  and  subtype  F.  tularensis  (12,  31)  and  was  used  as  the 
reference  methodology  to  evaluate  a  collection  of  88  wild-type 
and  4  well-characterized  F.  tularensis  isolates.  PFGE  was  per¬ 
formed  by  using  two  different  restriction  endonucleases  (Pmel 
and  BamHI).  In  these  PFGE  analyses,  F.  tularensis  strains 
SCHU  S4,  WY96,  and  LVS  were  used  as  subtype  A.I,  subtype 
A.II,  and  type  B  controls,  respectively.  Overall,  PFGE  analyses 
of  Pmel-digested  (Fig.  2)  and  BamHI-digested  (Fig.  3)  F.  tu¬ 
larensis  chromosomal  DNA  were  in  agreement,  identifying  56 
type  A.I,  15  type  A.II,  and  21  type  B  isolates.  PFGE  using 
Pmel  segregated  the  subtype  A.I  isolates  largely  into  two 
groups,  clades  A.I.a  and  A.I.b,  as  was  previously  defined  by 
Kugeler  and  associates  (18).  However,  PFGE  analysis  utilizing 
BamHI  did  not  group  the  subtype  A.I  isolates  into  the  same 
clade  A.I.a  and  A.I.b  subgroups.  A  detailed  analysis  of  the 
clustering  patterns  obtained  with  Pmel  and  BamHI  identified 
live  and  six  subpopulations  in  the  subtype  A.I  strains,  respec¬ 
tively.  The  PFGE  patterns  of  the  15  subtype  A.II  isolates  were 
more  divergent  than  those  of  the  subtype  A.I  strains,  as  pre¬ 
viously  reported  (18).  PFGE  analysis  with  both  Pmel  and 
BamHI  indicated  the  presence  of  six  discernible  subtype  A.II 
groups,  whereas  all  21  type  B  isolates  were  indistinguishable. 
For  F.  tularensis  subsp.  novicida  and  F.  philomiragia,  PFGE 
patterns  obtained  with  Pmel  and  BamHI  indicated  that  they 
were  both  outliers  and  distantly  related  to  F.  tularensis  type  A 
and  type  B  strains  (data  not  shown).  Together,  these  data 
validate  the  PFGE  patterns  obtained  with  the  use  of  cither 
Pmel  or  BamHI  to  identify  the  F.  tularensis  subspecies  or 
subtype. 

Subspecies  identification  and  subtyping  by  PCR.  The  CR10 
and  CR16  primers  were  utilized  in  a  PCR  to  type  a  collection 
of  wild-type  F.  tularensis  isolates.  Typing  results  obtained  with 
the  CR-derived  amplicons  were  compared  to  results  of  PFGE 
genotyping  with  Pmel  and  BamHI.  These  results  are  summa¬ 
rized  in  Table  2.  The  CR-derived  amplification  products  dif¬ 
ferentiated  the  F.  tularensis  strains  into  types  and  subtypes  that 
corroborated  the  results  of  genotyping  by  PFGE. 

The  majority  of  subtype  A.I  isolates  analyzed  with  primers 
CR10  C  and  CR10  S  generated  PCR  products  that  were  1.2  kb 
in  size  (43  strains).  An  amplicon  of  this  length  was  also  ob¬ 
served  for  the  subtype  A.I  prototype  strain,  SCHU  S4.  The 
remainder  of  the  subtype  A.I  isolates  (14  strains)  produced  a 
1.5-kb  amplicon,  indicating  that  this  locus  differed  by  approx¬ 
imately  300  bp  from  the  other  subtype  A.I  strains  (Table  2). 
Primer  pair  CR10  C  and  CR10  L  produced  a  1.3-kb  amplicon 
for  all  of  the  subtype  A.II  isolates  and  a  2.1-kb  amplicon  for  all 
of  the  type  B  strains  (Table  2).  For  F.  tularensis  subsp.  novicida , 
a  0.4-kb  PCR  product  was  produced  with  primers  CR10  C  and 
CR10  L,  and  no  amplicon  was  observed  with  primer  pair  CR10 
C  and  CR10  S. 

The  utilization  of  primers  CR16  C  and  CR16  S  produced 
either  a  2.3-kb  (54  strains)  or  a  2.0-kb  (3  strains)  amplicon  for 
the  F.  tularensis  subtype  A.I  isolates,  with  the  larger  PCR 
product  being  amplified  from  the  SCHU  S4  genome  (Table  2). 
For  the  subtype  A.II  strains,  primers  CR16  C  and  CR16  L 
produced  either  a  1.7-kb  or  a  1.4-kb  PCR  product,  with  the 
generation  of  the  larger  1.7-kb  amplicon  from  WY96  (Table 
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LVS 
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FTL_1 463 
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I 
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(2,333  bp) 
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1,799,532 
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FTT_1718c 

(ISFfuf) 


FTT_1719c 
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WY96 
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(ISFfuf) 
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\ 
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FTW  0083 
(fnuDI) 
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(781  bp) 


\ 
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I 
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FIG.  1.  Chromosomal  map  of  CR10  and  CR16  primer  pair  locations  in  F.  lularensis  subsp.  lularensis  subtypes  A.I  and  A.II  and  F.  tularensis 
subsp.  Iwlarctica  type  B.  The  chromosomal  regions  amplified  with  the  CR10  (A)  and  CR16  (B)  primer  sets  from  the  sequenced  reference  F. 
lularensis  strains  are  shown.  SCFTU  S4  (GenBank  accession  no.  AJ749949),  WY96-341S  (WY96)  (GenBank  accession  no.  CPOOfifiOS),  and  LVS 
(GenBank  accession  no.  AM233362)  represent  subtype  A.I,  subtype  A.II,  and  type  B  strains,  respectively.  The  resulting  antplicon  size  is  indicated 
below  the  appropriate  strain  in  parentheses.  The  forward  and  reverse  CR  (C,  S,  or  L)  primers  arc  denoted  by  black  arrows.  The  corresponding 
genes  in  the  various  strains  for  each  panel  arc  indicated  by  the  same  color  with  the  locus  identifier  number.  The  beginning  and  ending  nucleotide 
positions  for  the  mapped  loci  are  denoted  below  the  respective  chromosomal  regions.  Dashed  lines  designate  relative  locations  and  organizations 
of  the  corresponding  nucleotide  sequences  in  the  F.  tularensis  genomes. 


2).  All  type  B  isolates  produced  a  0.8-kb  amplicon  with  primers 
CR16  C  and  CR16  L  (Table  2).  No  amplicons  were  obtained 
with  either  primer  pair  CR16  C  and  CR16  S  or  CR16  C  and 
CR16  L  for  F.  tularensis  subsp.  novicida. 

Overall,  these  results  demonstrated  that  the  same  PCR 
product  profile  was  generated  with  all  four  CR  primer  pairs  for 


both  the  reference  and  wild-type  strains  of  the  same  subtype 
and/or  type.  To  reduce  the  number  of  PCR  samples  and  fur¬ 
ther  test  the  capability  of  the  DISA  assay,  the  CR10  and  CR16 
primer  pairs  were  reduced  to  only  two  sets  of  primers,  with 
each  set  containing  the  three  associated  primers  (CR10  C, 
CR10  S,  and  CR10  L  as  well  as  CR16  C,  CR16  S,  and  CR16  L). 
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Type 
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MO-B05 

2005 

AS 
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AJ 

AJ 
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A.I 

Al 
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A.I 
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Al 

AJ 
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A.I 

Al 
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A.I 

Al 
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2007 

Al 

AJ 

NE-1 02808 
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AJ 

Al 
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A.I 

Al 
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AJ 

AJ 
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AJ 

AJ 
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Al 
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Al 

AJ 
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AJ 

AJ 
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A.I 
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AS 

AJ 
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A.I 
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AJ 
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AJ 

Al 
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2005 

AJ 
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Al 

Al 
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A.I 
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AJ 
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Al 

AJ 
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AJ 

AJ 
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Al 

AJ 

NE-050707 

2007 

Al 

AJ 

NE-051807 
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Al 

Al 
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Al 

Al 

NE-073007 

2007 

AJ 

AJ 
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Al 

Al 
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2007 

Al 

AJ 
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Al 

Al 
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Al 

AJ 
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AJ 
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1999 

Al 

AJ 

WY-0705409 

2007 

AJ 

AJ 
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AJ 

Al 

NC-54558-01 

2001 

AJ 

Al 

NC-54S59-01 

2001 

AJ 

A.I 

UTPH-2173 

unknown 

AJ 

AJ 

NC-RabbilNC 

unknown 

AJ 

Al 

AK-1 100559 

2004 

AJ 

AJ 

NE-071605 

2005 

AJ 

A.I 

NE-062506 

2006 

AJ 

Al 

NE-0S2407 

2007 

AJ 

AJ 

NE-082007 

2007 

Al 

Al 

NE-050108 

2008 

AJ 

Al 

NE-052608 

2008 

AJ 

AJ 

NE-062008 

2008 

AJ 

Al 

NE-071408 

2008 

Al 

AJ 

NE-110607 

2007 

AJ 

Al 

NE -061 608 

2008 

Al 

AJ 

FIG.  2.  Dendrogram  of  PFGE  patterns  obtained  with  Pmel-digested  F.  lularensis  subsp.  lularensis  subtype  A.I  and  the  apparent  clade  A.I.a  and 
A.I,b  subpopulations  (A)  as  well  as  F.  lularensis  subsp.  tularensis  subtype  A.II  and  F.  lularensis  subsp.  holarclica  type  B  (B).  A  representative 
Pmel-digested  F.  tularensis  A.I  strain  from  each  of  the  five  clusters  is  shown  in  the  upper  part  of  B  for  comparison.  Distance  is  displayed  in  relative 
units. 
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Year  PFGE  DISA 


Strain 

Infected 

Type 

Type 

NE-111407 

2007 

At 

Al 

NC-48620-97 

1997 

A.l 

Al 

SCHU  S4 

1941 

A.I 

A.I 

AK-1 133496 

2003 

Al 

Al 

NE-1 10607 

2007 

Al 

Al 

WY-05F1990 

2005 

All 

All 

WY-05W9954 

2005 

All 

All 

WY-06W9410 

2006 

All 

All 

WY-07F13554 

2007 

All 

All 

UT-80402860 

2004 

All 

All 

WY-BT254 

unknown 

All 

All 

ATCC  6223 

unknown 

All 

All 

WY-03W10146 

2003 

All 

All 

UT-70102163 

2001 

All 

All 

WY-HumanNC 

unknown 

All 

All 

WY-BT324 

1996 

All 

All 

WY96-3418 

1996 

All 

All 

WY-06F12590 

2006 

All 

All 

WY-00W41 14 

2000 

All 

All 

WY-BT188 

unknown 

All 

All 

LVS  (ATCC  29684) 

unknown 

B 

B 

WY-961 94280 

1996 

B 

B 

WY-9868529 

1998 

B 

B 

WY-OvineNC 

unknown 

B 

B 

MO-MS304 

unknown 

B 

B 

MO-HumanNC 

unknown 

B 

B 

FR-SS 

unknown 

B 

B 

FR-LR 

1993 

B 

B 

NE-061705 

2005 

B 

B 

MO-E05 

2005 

B 

B 

MO-G05 

2005 

B 

B 

UT-70001092 

2000 

B 

B 

UT-70101201 

2001 

B 

B 

WY-BT235 

unknown 

B 

B 

WY-BT323 

1996 

B 

B 

NE-072606 

2006 

B 

B 

WY-05W8084 

2005 

B 

B 

WY-05W13621 

2005 

B 

B 

WY-05W13374 

2005 

B 

B 

NE-1 11 108 

2008 

B 

B 

NE-073009 

2009 

B 

B 

These  tests  produced  the  same  PCR  product  profiles  as  those 
of  the  four  primer  pairs  hut  with  half  as  many  reactions. 

Amplification  specificity  and  sensitivity  of  CR  primers.  The 
CR  primer  specificity  for  F.  tularensis  was  evaluated  by  using 
an  exclusivity  panel  (n  =  24)  (Table  1)  that  included  chro¬ 
mosomal  DNA  from  the  near  neighbors  F.  philomiragia  and 
Wolbachia  persica.  Neither  the  CR10  nor  CR16  primer  pairs 
resulted  in  any  prominent  PCR  products  for  these  bacterial 


samples,  with  the  exception  of  W.  persica,  in  which  a  250-bp 
amplicon  was  generated  by  use  of  primers  CR10  C  and 
CR10  L. 

The  applicability  of  PCR  amplification  with  the  CR  primers 
in  environmental  studies  to  detect  the  presence  of  F.  tularensis 
was  examined  by  using  arthropods,  pond  water,  and  soil.  Fifty- 
four  ticks,  including  Demiacentor  andersoni,  Dermacentor 
variabitis,  and  Amhlyomma  americanum ,  species  known  to 
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Type 
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Al 
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Al 
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Al 
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Al 

Al 
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Al 

Al 
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A.I 

Al 

NE-052608 
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Al 

Al 
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A.I 

A.I 
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2008 

Al 

A.I 

NE-071408 

2008 

Al 

Al 

WY-07F11688 
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Al 

A.I 

WY-07F10961 

2007 

A.I 

A.I 

NE-110607 

2007 

A.I 

A.I 

WY-07F8058 

2007 

A.I 

A.I 

WY-07F109C8 

2007 

Al 

Al 

MO-MS1349 

unknown 

Al 

Al 

NE-100903 

2003 

Al 

A.I 

NE-07190S 

2005 

A.I 

Al 

NE-121905 

2005 

Al 

Al 

MO-B05 

2005 

Al 

Al 

MO-C05 

2005 

Al 

Al 

NE-050206 

2006 

A.I 

Al 

NE-073007 

2007 

Al 

Al 

NE-092107 

2007 

A.I 

Al 

NE-1 20707 

2007 

Al 

Al 

WY-06F14621 
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Al 

Al 

WY-0705409 
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Al 

Al 

NE-1 11407 

2007 

A.I 

Al 

NE-102908 

2008 

A.I 

Al 

FIG.  3.  Dendrogram  of  PFGE  patterns  obtained  with  BamHI-digested  F.  tularcnsis  subsp.  tularensis  subtype  A.I  (A)  as  well  as  F.  tularensis 
subsp.  tularensis  subtype  A.II  and  F.  tularensis  subsp.  holarctica  type  B  (B).  A  representative  BamHI-digested  F.  tularensis  A.I  strain  from  each  of 
the  sue  dusters  is  shown  in  the  upper  part  of  B  for  comparison.  Distance  is  displayed  in  relative  units. 
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Strain 

Year 

Infected 

PFGE 

Type 

DISA 

Type 

SCHU  S4 

1941 

Al 

At 

NE-1 10607 

2007 

A.I 

A.I 

AK-1 100558 

2004 

A.I 

At 

WY-07F10908 

2007 

A.I 

A.I 

NE-092107 

2007 

Al 

A.I 

NE-1 11 407 

2007 

Al 

Al 

UT-70102163 

2001 

All 

All 

WY-07F 13554 

2007 

All 

All 

WY-05W9954 

2005 

All 

All 

WY-06W9410 

2006 

Alt 

All 

UT-80402860 

2004 

All 

All 

WY-BT254 

unknown 

All 

All 

ATCC  6223 

unknown 

All 

All 

WY-03W10146 

2003 

All 

All 

WY-00W4114 

2000 

All 

All 

WY-BT188 

unknown 

All 

All 

WY-HumanNC 

unknown 

All 

All 

WY-BT324 

1996 

All 

All 

WY96-3418 

1996 

All 

All 

WY-05F1990 

2005 

All 

All 

WY-06F12590 

2006 

All 

All 

LVS  (ATCC  29884) 

unknown 

B 

B 

WY-96 194280 

1996 

B 

B 

WY-9868529 

1998 

B 

B 

WY-OvineNC 

unknown 

B 

B 

MO-MS304 

unknown 

B 

B 

MO-HumanNC 

unknown 

B 

B 

FR-SS 

unknown 

B 

B 

FR-LR 

1993 

B 

B 

NE-061 705 

2005 

B 

B 

MO-E05 

2005 

B 

B 

MO-G05 

2005 

B 

B 

UT-70001092 

2000 

B 

B 

UT-70101201 

2001 

B 

B 

WY-BT235 

unknown 

B 

B 

WY-BT323 

1996 

B 

B 

NE-072606 

2006 

B 

B 

WY-05W8084 

2005 

B 

B 

WY-05W13621 

2005 

B 

B 

WY-05W13374 

2005 

B 

B 

NE-111108 

2008 

B 

B 

NE-073009 

2009 

B 

B 

FIG.  3 — Continued. 


transmit  F.  tularensis  in  North  America  (26),  were  assessed  for 
the  presence  of  F.  tularensis. 

PCR  with  the  CR  primers  revealed  that  one  D.  andersoni 
tick  and  one  D.  variabilis  tick  produced  products  consistent 


with  the  presence  of  F.  tularensis,  whereas  the  remainder  of  the 
ticks  did  not  produce  any  amplicons.  Importantly,  the  PCR 
product  sizes  obtained  with  the  CR  primers  provided  subspe¬ 
cies  and  subtype  details.  Primer  pair  CR10  C  and  CR10  L  and 
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TABLE  2.  Comparison  of  typing  results  obtained  by  PCR  with  the  CR  primers  and  by  PFGE  with  Pmel  or  BamHI  for  F.  tularensis  subsp. 
lularensis  (subtypes  A.I  and  A.I1),  F.  tularensis  subsp.  Iiolarciica  (type  B),  and  F.  lularensis  subsp.  novicida  isolates 


No.  of  F.  lularensis 
isolates 

CR  primer  pair 

Amplicon  size  (kb) 

Type  determined  by: 

PCR-based  typing 
with  CR  primers 

PFGE-based  typing 
with  Pmel 

PFGE-based  typing 
with  BamHI 

43° 

CR10  C  +  CR10  S 

1.2 

A.I 

A.I 

A.I 

14 

CR10  C  +  CR10  S 

1.5 

A.I 

A.I 

A.I 

54“ 

CR16  C  +  CR16  S 

2.3 

A.I 

A.I 

A.I 

3 

CR16  C  +  CR16  S 

2.0 

A.I 

A.I 

A.I 

15* 

CR10  C  +  CR10  L 

1.3 

A.II 

A.II 

A.II 

14" 

CR16  C  +  CR16  L 

1.7 

All 

A.II 

A.II 

1 

CR16  C  +  CR16  L 

1.4 

A.II 

A.II 

A.II 

21c 

CR10  C  +  CR10  L 

2.1 

B 

B 

B 

21“ 

CR16  C  +  CR16  L 

0.8 

B 

B 

B 

1" 

CR10  C  +  CR10  L 

0.4 

F.  tularensis  subsp. 

F.  lularensis  subsp. 

F.  lularensis  subsp. 

novicida 

novicida 

novicida 

1" 

CR10C  +  CR10S 

1" 

CR16  C  +  CR16  L 

1" 

CR16  C  +  CR16  S 

“  Includes  F.  lularensis  subsp.  lularensis  strain  SCHU  S4  (subtype  A.I). 

*  Includes  F.  lularensis  subsp.  lulurensis  strain  WY96-3418  (subtype  A.II), 
'  Includes  F.  lularensis  subsp.  Iiolarciica  LVS  (type  B). 

J  F.  lularensis  subsp.  novicida  Utah  112  (ATCC  15482). 


primer  pair  CR16  C  and  CR16  L  produced  amplicons  indicat¬ 
ing  the  presence  of  F.  lularensis  subsp.  tularensis  subtype  A.II 
in  the  D.  andersoni  tick.  Amplification  products  derived  with 
primer  pair  CR10  C  and  CR10  L  and  primer  pair  CR16  C  and 
CR16  L  and  DNA  isolated  from  the  D.  variabilis  tick  were 
characteristic  of  F.  lularensis  subsp.  holarctica  (type  B).  The 
amplicons  derived  from  these  two  ticks  were  sequenced,  and 
the  data  corresponded  to  the  appropriate  F.  tularensis  subtype 
and/or  type  as  determined  by  the  DISA  assay.  Results  from 
two  commercially  available  PCR-based  test  kits  (Idaho  Tech¬ 
nology,  Inc.,  Salt  Lake  City,  UT)  that  detect  F.  tularensis  were 
also  positive  using  DNA  extracts  from  these  two  ticks.  How¬ 
ever,  the  commercial  assays  could  not  type  the  F.  tularensis 
isolates.  No  viable  F.  tularensis  was  recovered  from  these  two 
tick  specimens  by  culturing  on  chocolate  agar  plates,  and  no 
other  environmental  samples  were  positive  with  either  the 
commercial  or  CR  primer  pairs.  Together,  these  data  demon¬ 
strated  that  the  CR  primers  were  specific  for  F.  tularensis  and 
capable  of  detection  and  subspecies  or  subtype  identification 
for  both  the  clinical  and  environmental  isolates  examined. 

As  a  control  to  verify  the  presence  of  Francisella-Yike  endo- 
symbionts  in  the  ticks,  PCR  amplifications  were  performed  by 
using  primers  known  to  amplify  target  genes  in  Francisella- like 
bacteria  (3,  11,  24,  30).  The  Francisella-Uke  endosymbiont  W. 
persica  served  as  a  positive  control.  In  addition,  commercially 
available  test  kits  (Idaho  Technology,  Inc.)  that  detect  only  F. 
tularensis  subspecies  were  utilized  in  these  assessments.  Ticks 
that  produced  PCR  products  with  the  primers  designed  to 
amplify  Francisella-Uke  bacteria  were  further  evaluated  with 
the  commercial  F.  tularensis-specific  assays  as  well  as  with  the 
DISA  assay.  If  no  amplicons  were  produced  with  the  F.  tular- 
e/u/r-specific  assays  but  PCR  products  were  produced  with  the 
primers  for  Francisella-Uke  bacteria,  the  respective  amplicons 
were  sequenced.  The  subsequent  nucleotide  sequences  were 
compared  to  those  deposited  in  the  databases,  confirming  the 
presence  of  Francisella-Uke  endosymbionts  and  the  absence  of 
F.  tularensis  or  F.  philomiragia  in  the  associated  ticks.  These 


findings  demonstrated  that  the  CR  primers  did  not  amplify 
DNA  from  Francisella-Uke  endosymbionts  known  to  exist  in 
ticks. 

The  limit  of  detection  for  the  PCR  amplification  of  F.  tula¬ 
rensis  with  the  CR  primers  under  standard  PCR  conditions 
described  in  Materials  and  Methods  was  determined  to  be 
approximately  20  copies  of  the  1.9-Mbp  F.  tularensis  chromo¬ 
some.  Therefore,  PCR  with  the  CR  primers  provided  a  detec¬ 
tion  of  F.  tularensis  at  a  level  naturally  present  in  environmen¬ 
tal  sources  and  achieved  an  identification  of  the  subspecies  or 
subtype. 

Sequence  analysis  of  amplicons.  Representative  PCR  prod¬ 
ucts  obtained  with  the  CR  primers  for  each  size  group  were 
cloned  and  sequenced  to  determine  content.  The  DNA  se¬ 
quences  of  the  PCR  products  generated  with  the  CR10  and 
CR16  primers  for  these  strains  were  in  agreement  with  the 
expected  lengths  (Table  2)  and  mapped  loci  (Fig.  1).  The  size 
difference  observed  for  the  1.2-kb  and  1.5-kb  amplicons  in  the 
F.  tularensis  subtype  A.I  clade  obtained  with  primers  CR10  C 
and  CR10  S  was  due  to  a  312-bp  deletion  in  the  smaller  PCR 
product  (Table  3). 

With  primer  pair  CR16  C  and  CR16  S,  either  a  2.3-kb  or  a 
2.0-kb  amplicon  was  generated  for  the  subtype  A.I  strains, 
while  primers  CR16  C  and  CR16  L  resulted  in  either  a  1.7-kb 
or  a  1.4-kb  PCR  product  for  the  subtype  A.II  isolates.  In  both 
cases,  the  smaller  amplification  product  was  again  due  to  a 
312-bp  deletion  (Table  4).  The  312-bp  deleted  region  in  the 
smaller  CR10-  and  CR16-derived  amplicons  correlated  with 
the  middle  region  of  the  IS Flu  I  transposase  gene.  Single-nu¬ 
cleotide  polymorphisms  (SNPs)  were  also  noted  for  some  of 
the  CR-derived  PCR  products,  as  indicated  in  Tables  3  and  4. 

DISCUSSION 

As  with  other  virulent  pathogens,  such  as  Bacillus  anthra- 
cis  and  Yersinia  pestis,  F.  tularensis  is  highly  monomorphic, 
with  greater  than  99%  sequence  similarity  between  type  A 
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TABLE  3.  Sequence  analysis  of  representative  amplicons  obtained  by  PCR  with  the  CR10  primer  set  relative  to 

F.  tularensis  SCHU  S4,  WY96-3418,  or  LVS” 


F.  tularensis  isolate 

Type 

CR10  primer  pair 

Amplicon  size  (kb) 

Insertion 

SNP6 

SCHUS4 

A.I 

C  +  S 

1.2 

NE-061598 

A.I 

c  +  s 

1.2 

NC-52797-99 

AI 

c  +  s 

1.2 

nt  896/A— T;  nt  9557T— C 

NE-071905 

AI 

c  +  s 

1.5 

312  bp  at  nt  539-850” 

nt  512/T-C 

WY-0705409 

AI 

c  +  s 

1.5 

312  bp  at  nt  539-850” 

nt  1 135/A— »G;  nt  1244/T— C 

WY96-3418 

All 

C  +  L 

1.3 

WY-0OW4114 

A.II 

C  +  L 

1.3 

nt  755/T— C 

WY-HumanNC 

All 

C  +  L 

1.3 

LVS 

B 

C  +  L 

2.1 

“  Nucleotide  (nt)  sequences  for  reference  strain  SCHU  S4,  WY96-3418,  and  LVS  amplicons  were  obtained  from  GenBank  accession  no.  AJ749949,  CP000608,  and 
AM233362,  respectively. 

'■  SNP,  single-nucleotide  polymorphism. 

”  Derivative  of  IS  Flu  I. 


and  type  B  strains  (27, 28).  Genetically  monomorphic  patho¬ 
gens  may  have  recently  passed  through  a  population  bottle¬ 
neck,  reducing  or  abolishing  genetic  diversity  (1).  Despite 
the  considerable  conservation  of  overall  genomic  content 
among  F.  tularensis  subspecies,  the  organizations  of  discrete 
genetic  elements  differ  (8).  The  observed  translocations  and 
rearrangements  appear  to  be  mediated  by  the  movement  of 
IS  elements,  most  notably  by  the  highly  abundant  IS Ftul 
gene  (27).  Seven  different  types  of  IS  elements  have  been 
found  within  SCHU  S4  and  the  closely  related  subtype  A.I 
strain  NE-061598,  including  IS  Ftul  (50  copies);  IS  Ftu2  (16 
copies);  lSFtu3  and  IS Ftu6  (3  copies  each);  and  IS Ftu4, 
ISFtuS,  and  lSSodl3  (1  copy  each)  (23).  IS  elements  appear 
to  be  the  primary  driving  force  for  the  evolution  of  the 
Francisella  chromosome,  contributing  to  genome  decay  and 
subsequent  intracellular  dependence. 

Although  there  are  similar  numbers  of  IS  Ftul  genes  in  F. 
tularensis  subsp.  tularensis  (type  A)  and  F.  tularensis  subsp. 
holarctica  (type  B)  and  many  more  copies  of  IS Ftu2  in  type 
B  than  in  type  A  strains,  genomic  rearrangements  in  type  B 
strains  are  much  less  frequent  (32).  The  low  level  of  genetic 
diversity  and  high  frequency  of  transposable  elements  in  the 
type  B  clade  support  the  hypothesis  that  F.  tularensis  subsp. 


holarctica  emerged  after  F.  tularensis  subsp.  tularensis  and, 
therefore,  is  more  clonal.  In  addition,  the  wider  geographic 
dispersion  of  F.  tularensis  subsp.  holarctica  supports  the 
correlation  between  this  pathogen  either  having  a  higher 
capacity  to  survive  in  different  environmental  niches  than  F. 
tularensis  subsp.  tularensis  or  causing  less  morbidity  in  the 
associated  host.  The  availability  of  a  rapid  PCR-based 
method  for  examining  additional  wild-type  isolates  may  pro¬ 
vide  new  insights  into  the  mechanism  for  the  observed  geo¬ 
graphic  associations  between  these  two  subspecies. 

Despite  the  subtype  A.I  and  A.II  strains  having  similar 
numbers  of  IS  elements  (6,  20),  there  is  considerable  overall 
variation  in  chromosome  architectures  (18,  23).  The  current 
study  revealed  that  a  312-bp  insertion  or  deletion  derived 
from  IS  Ftul  was  responsible  for  the  size  heterogeneity 
within  the  CR-derived  amplicons  in  the  subtype  A.I  and  A.II 
strains  but  that  discrimination  between  these  two  clades  was 
still  apparent.  In  Mycobacterium,  Salmonella,  and  other 
pathogens  the  locations  of  IS  elements  have  been  sufficiently 
stable  for  use  as  markers  in  restriction  fragment  length 
polymorphism  studies  and  typing  for  epidemiological  pur¬ 
poses  (21).  In  Francisella,  the  spatial  location  of  IS  elements 
and  their  involvement  in  the  plasticity  of  the  genome  can 


TABLE  4.  Sequence  analysis  of  representative  amplicons  obtained  by  PCR  with  the  CR16  primer  set  relative  to 

F.  tularensis  SCHU  S4,  WY96-3418,  or  LVS” 


F.  tularensis 
isolate 

Type 

CR16 
primer  pair 

Amplicon 
size  (kb) 

Insertion 

Deletion 

SNP” 

SCHU  S4 

A.I 

c  +  s 

2.3 

NE-061598 

A.I 

c  +  s 

2.3 

WY-07F10908 

A.I 

c  +  s 

2.0 

T  at  nt  1193 

312  bp  at  nt  1282-15936 

nt  1202/T— C 

NE-052608 

A.I 

c  +  s 

2.0 

312  bp  at  nt  1282-15936 

WY-07F8058 

AI 

c  +  s 

2.0 

T  at  nt  1193 

312  bp  at  nt  1282-15936 

nt  1202/T— C 

WY96-3418 

A.II 

C  +  L 

1.7 

WY-OOW4114 

All 

C  +  L 

1.7 

WY-HumanNC 

All 

C  +  L 

1.7 

WY-07F13554 

All 

C  +  L 

1.4 

312  bp  at  nt  691-1002* 

nt  567/C— T;  nt  643/A— G;  nt  1126/G— C;  nt 

1228/T— C;  nt  1389/T— C;  nt  1650/T— C 

LVS 

B 

C  +  L 

0.8 

WY-05W8084 

B 

C  +  L 

0.8 

UT-70101201 

B 

C  +  L 

0.8 

nt  290/A— T 

*  Nucleotide  sequences  for  reference  strain  SCHU  S4,  WY96-3418,  and  LVS  amplicons  were  obtained  from  GenBank  accession  no.  AJ749949,  CP000608,  and 
AM233362,  respectively. 

'’Derivative  of  IS  Ftul. 

‘  SNP,  single-nucleotide  polymorphism. 
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also  be  exploited.  Our  future  studies  will  determine  whether 
the  DISA  method  will  provide  a  means  to  detect  and  mon¬ 
itor  the  emergence  of  new  substrains  of  this  highly  infectious 
pathogen. 

Type  A.I,  A.II,  and  B  strains  of  F.  tularensis  are  present  in 
North  America,  with  hard  ticks  being  considered  the  major 
biological  vector  of  tularemia  in  this  region  (26).  Since  F. 
tularensis  survival  is  sustained  throughout  the  tick’s  life  cycle 
from  larva  to  nymph  to  adult,  and  since  tick  morphogenesis 
requires  a  blood  meal,  pathogen  transmission  to  a  new  host 
may  occur  at  each  stage.  Therefore,  a  single  infected  tick 
may  be  responsible  for  establishing  a  new  or  recurring  focus 
of  disease  in  the  environment.  The  multidirectional  nature 
of  zoonosis,  the  presence  of  Francisella-like  endosymbionts 
in  ticks,  and  the  broad  host  range  of  F.  tularensis  have 
complicated  epidemiological  investigations  into  diseases 
caused  by  this  pathogen.  A  field-deployable  assay  would  be 
useful  for  the  characterization  of  infection  outbreaks  and 
the  identification  of  their  point  source  or  location.  The 
PCR-based  DISA  method  targeting  IS  elements  provides  a 
reliable  approach  for  the  rapid  evaluation  of  ticks,  animals, 
and  humans  potentially  infected  with  F.  tularensis.  Further¬ 
more,  an  accurate  surveillance  of  F.  tularensis  vectors  and 
reservoirs  will  provide  a  better  understanding  of  this  patho¬ 
gen’s  biology. 

Molecular-based  assays  for  the  detection  of  infectious 
diseases  have  continued  to  evolve  in  order  to  address  the 
problem  of  false-positive  results  due  to  the  cross-hybridiza¬ 
tion  of  primers  with  unrelated  sequences  with  unexpected 
homology.  Some  of  the  most  effective  changes  have  been  the 
incorporation  of  multiple  targets  into  the  assay,  with  the 
criterion  for  positivity  being  the  amplification  of  more  than 
one  target.  Although  several  real-time  PCR  assays  have 
been  developed  to  differentiate  isolates  of  F.  tularensis  that 
arc  of  types  A  and  B  (19)  or  types  A.I  and  A.II  (22),  a  more 
comprehensive  method  was  needed  to  simultaneously  and 
accurately  genotype  this  species.  The  DISA  assay  offers  this 
capability  and  requires  only  two  sets  of  primers,  unlike  typ¬ 
ing  by  PCR-based  multilocus  variable-number  tandem  re¬ 
peat  analysis.  Furthermore,  the  DISA  method  takes  advan¬ 
tage  of  the  presence  of  multiple  inherent  targets  within 
Francisella,  and  this  approach  will  be  used  for  the  future 
development  of  a  real-time  molecular  assay. 
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